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ABSTRACT 


This report describes an unsteady pressure dleteibution cus 
to a blaet wave ciffracting eround a SPARTAN missile ascembly. 
The pressure dava were obtained from a series of five teets par- 
formed in the DASACON Conical Shock Tube Facility iocated at tse 
Naval Weapons Laboratory, Dahlgren, Virginia. These tests were 
conducted during the period 17 April 1972 to 8 May 1972. During 
these tests the aissile assembly wus subjected to incident blast 
waves which had peak overpressures of from 2.9 pei to 11.8 pei 
and corresponding positive overpreseurs duretiou oc? from 
380 milliseconds to 444 milliseconds. 


The report describes the pressure data for each test and 
the empirical function used to represent thes data. It then 
describes the method of integrating the empirical function et 
given times for the missile assembly sections of interest. Tue 
results of these calculations for all five tests exe given at 
selected times. The calculation period covers arproximately 
seven milliseconds teginning at the time the blaet wave iret 
encounters the missile assembly. These resuits ere given e5 
force ys time plots ita Appendix D. 


FOREWORD 


This work was performed by the Naval Weapons Laboratory (NWL), 
Dahlgren, Virginia, for McDonnell-Douglas Astronautics Company (MDAC) 
in accordance with MDAC test control drawing number 1T4-7031, 

29 March 1972. Funding was provided by the U. S. Army Safeguard 
Systems Commander under MIPR No. A31699-23-V180, PRON No. OR-23-V180, 
and CMC Code 527A.000. 


This report was reviewed by J. J. Yagla and F. F. Churchill of 
the Test and Evaluation Department. 


Released by: 


ELE Clengbery 


L. A. CLAYBERG, Head 
Test and Evaluation Department 


ii 


FOREWORD 


ABSTRACT . 
I. INTRODUCTION ¢ dag , ; 5 
II. DESCRIPTION OF PRESSURE DISTRIBUTION Fels. Hel as 


A 
B. 
C. 


IV. DESCRIPTION OF FORCE CALCULATION 


a 
B. 
Cc. 


III. DESCRIPTION OF PRESSURE FIT ..... 


CONTENTS 


eo. e oe «8 e e o @ . ee e e @ . @ 


SHOCK DIFFRACTION PROCESS .... 
MSASUREMENTS ON FINLESS SIDE . 
MEASUREMENTS ON FIN SIDE... . 


CYLINDRICAL SECTIONS. ....- 
FRUSTRUM SECTIONS ee @ . o ee e 


DSSCRIFTION OF COMPUTER PROGRAM We. des oe 


Mt ° RESULTS e 8 . eo sere <8 8. e oo ee e@# © @ e@ e eo 8 


VI. CONCLUSIONS AND RECOMMENDATIONS ........-s 


REFERENCES ... . fe a Rea se oP je) Bhp Sis “ere ces 


FIGURES: 


i. 
2. 


SPARTAN Missile Assembly 

Conical Shock Tube Facility 

Suspended Missile Assembly 

Pressure Transcucer Locations on Surface of 
Missile Assembly, Tests 114-115 

Pressure Transducer Locations on Surface of 
Missile Assembly, Tests 116-117-118 

Shock Wave Diffraction Around a Cylinder 
Coordinates on Surface of Missile Assembly 


F 


ron 
HO ONT 


ral oe 


SPS RS tet 
OFAANO DW 


SPARTAN 


SPARTAN 


Effect of Pressure Fit When Shock is Between Two 


Transducer Locations 

Intersection of Shock with Missile Assembly 
Iutegration Area for a Cylindrical Section 
actegration Areas of a Frustrum Section for 
Stages of Engulfment 

Flew Chart for Computer Program, PREDIN 
Pressure vs Angle on Control Section - Time 
Pressure vs Angle on Control Section - Time 
Pressure vs Angle of Control Section - Time 
Pressure vs Angle of Control Section - Time 
Force Ratio vs Time for Control Section 


Sections 


Ver ious 


Unit Normal on the Surface of a Fustrum Section 


Summary of Test Conditions 


Values of Summation Limits for Various Calculation 


Times 
1ii 


CONTENTS (Continued) 


APPENDICES 
A. Smoothed Pressure-Time Histories 
B. Pressure ve Time ca Fin Side of Control Section 
C. Pressure Function Coefficients 
D. Comparisons o? Curve Fits with Experimental Data 
E. Force-Time Histories 
F. Distribution 


iv 


I. INTRODUCTION 


The analysis described in this report was performed for 
McPonnell-Douglas Astronautics Company in connection with a series 
of five air blast loading tests on a SPARTAN missile assembly. 
This assembly, shown in Figure 1, consisted of a nose fuiring, the 
control section, the guidance section, and a dummy warhead. The 
shaded sections are the ones on which force calculations were 
performed. 


The missile assembly was tested in the 22 foot diameter test 
area (test aren 3) of the conical shock tube (DASACON) at the 
Naval Weapons [aboratory (see Figure 2). These tests are desig- 
nated by, DASACON test numbers 114-118, and are described in refer- 
ence (a)*. As shown in Figure 3, the missile assembly was suspended 
from the top of the test area in a nose down position. ‘The peak 
overpressures and the positive durations of the incident blast waves, 
along with the amblent temperature and pressure, are listed for 
ench test in Table 1. 


The objectives of these tests were to measure accelerations 
on portions of the assembly's primary structure, and overpressures 
on the assembly’ surface. This report is concerned with the 
surface press e di: ‘r‘bution. The pressure data were fitted to 
the eqration, 


K OL t-1 
p(s,*) = Y Y by ys cos [(j-1)¢ ] (1) 
$e). 4%: x 


p is pressure, 


ogitudinel distance along the surfacé measured from 


is 
the assembly tip, 


8 is circumferential angle measured from the forward stagna- 
tion line on assembly surface, 


bij are coefficients determined by fitting the equation to the 
experimental data. 


The data fitting was performed at given times for a period of 
approximately ceven milliseconds. By means of equation i; the 
pressure can be determined at any point on the surface of the 
major sections of the missile assembly for the given times. 


o_ = oe B= ew em mw mmm 


line references are located on page 38. 
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Such an expression is especially useful in calculating pres- 
sure forces ou the assembly sectiozs. A particular application for 
the DASACON tests was to integrate the pressure over the three 
major sections of the missile assembly. These integrations will 
be called "forceg" throughout this report although it is not quite 
proper to do so.“ Thus, the word "force", without any qualifica- 
tion, is to be taken as the integral of the pressure. Force vs 
time data was used by MNAC as a characteristic of the structural 
loading of the missile assembly produced by the DASACON environ-~ 
mert. ‘These data were compared with similar data produced by 
other environments. 


The objectives of this report are to deecribe the pressure 
data taken during the DASACON tests ang to show how these data 
were used to obtain the force calculations. The pressure deta are 
described in Section IZ and the method of fitting the data to 
equation (i) is discussed in Section III. The integration of 
equation (1) is discussed in Section IV, where the following 
assumptione are made: 


a. The pressure dietribution is symmetrical, i.e., the 
effects of the fin located on the control section (See Figure i) 
are ignored. 


b. The shock wave is plane as it diffcactsa around a giver 
section. 


Section V discusses the results of the force calculations and 
estimates the maximum effects of including measurements on the 
fin side in the analysis. Tables of the tj; coefficients in 
equation (1) are given in Appendix C for each calculation time 
of each test. Comparisons of the results of the pressure fits 
and experimental data are given in Appendix D. Section VI gives 
conclusions and recommendations. 


II. DESCRIFTION OF PRESSURE DISTRIBUTION 


Figures 4 and 5 show the locations of the pressure trans- 
ducers om the migsile assembly surface during the DASACON tests. 
The differences between these figures are due to a 180° rotation 
ef the assembly abcut its axis after the second test of the 


“me integral of the pressure over the assembly surface has the 
dimensions of a force, wut does not specify any direction. Thus 
it cannot be properly regarded as a vector quantity such as 
force. . ‘ sient 
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series (test No. 115). The number of pressure transducers on the 
surface of the assembly was limited to 30. Therefore, in order to 
define the pressure distribution more accurately, most of these 
gauges were placed on one side of the assembly. ‘Transducers were 
placed on the side opposite the fin at intervals of 30° around the 
circumferences of the sections of interest (see Figure 1). Three 
transducers were placed at stations 1 and , so that data from the 
boundaries of the integration region could be obtained. Three 
transducers, designated by the letter F were placed at station 3 
(control section) on the fin side to determine the maximum effects 
of the fin. 


The data from these measurements were recorded on magnetic 
tape (reference (a)). In order to have smooth data for the calcula- 
tions, the tapes were played back through low pass, linear phase 
analog filters with a cutoff frequency of 8 KHz. The data were 
then digitized, and plotted by a CALCOMP plotter. The resulting 
picts were carefully examined and then smoothed by hand. The 
smoothed data were digitized on punched cards at a rate of 
1000 samples per second for use in a computer. 


The smoothed data from each test are given in Apperdix A for 
the finless side of the missile assembly. Unsmoothedi data for 
the fin side are given in Appendix B. The appendices contain 
graphs of pressure vs time at the given circumferential angles. 
The station number for each curve is listed on the right side of 
each figure. The zero time in the figures of Appendix A is the 
time that the shock first touches the missile assembly. The 
symbol, T,, located at the lower left of the figures ir Appendix B, 
designates an arbitrary common time. 


A. Shock Diffraction Process 


Before the pressure records are discussed, a brief 
description of the shock diffraction process will be givea. 
Since the cross section at each axial position of the missile 
assembly is circular, the shock diffraction at a given axial 
position is similar to the diffraction of a shock around a 
cylinder. 


Shadowgrams of the diffraction of a shock wave around a 
cylinder are given in reference (b). Figure 6, which shows the 
diffraction at various stages of engulfment, is based on these 
shadowgrams. Figure 6a shows a regular reflection of the shock 
wave on the front of the cylinder surface. Regular reflection 
occurs until the angle, a, between the plane of the shock and the 
tangent to the surface is some critical value Qo, (reference (c)). 
This angle, Q.,, is dependent upon the incident shock strength. 
After the shock has reached the critical angle, a Mach stem is 
formed (Figure 6b). 
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This Mach stem weakens as it diffracts around the rear of 
the cylinder (Figure 6c). Mach stems for each side of the cylinder 
intersect and cross each other at @ = 180°. After crossing each 
other, they propagate back toward the front of the cylinder 
(Figure 64). As the shocks continue to propagate toward the front 
ier, the flow behind them is characterized by separation 
- formation. The separation is due to the adverse pressure 
gradient, caused by the shocks, on the beundary layer flow 
(Figure 6e). 


B. Measurements on Finless Side 


Figures A-8 thru A-10 of Appendix A sre typical of records 
from the front of the missile assembly, where regular reflection 
and initial formation of the Mach stem takes place. These records 
show an instantaneous rise to a peak followed by a rapid decay and 
then by a much more gradual decay. Figures A-1l through A-14 are 
typical of records from the rear of the assembly. The curves of 
Figures A-ll thru A-13 ahow an instantaneous rise to a peak followed 
by a decay to a minimum value and then a second rise. The time 
between the pressure rises decreases as the circumferential angle 
increases. The second pressure rise is caused by the shock from 
the other side of the cylinder being propagated toward the front. 
Figure A-14 shows the pressure records from the 9 = 180° locations. 
These curves show peaks which heave values high=r than those of the 
90° < 9 <150° locations. The higher pressures are caused by the 
shock on one side of the cylinder intersecting tae shock from ths 
other side. 


The curves of Appendix A for a given 9 location show that 
the pressure variation with the longitudinal distance on the 
frustrum sections (sections, one, two, and three shown in Figure 1) 
is quite pronounced on the rear of the missile assembly. These 
variations are lergely caused by the fact that the shock reaches 
the smaller diameter sections at later times then the larger ones 
and by the shork's engulfing the smaller diameter sections in less 
time. Thus, the pressure maxims and minimums occur at different 
times for the different frustum sections. 


C. Measurements on Fin Side 


Tests 115 and 116 were conducted at incident peak shock 
overpressures of approximately 12 psi. Since the missile assembly 
was rotated 180° after test 115, pressure measurements at locetions 
of 30°, 60°, 75°, 105°, 120° and 150° were obtained on the fin side 
of station 3 for essentially the same incident pressure. Pressure- 
time plots from these measurements are shown in Appendix B in order 
of increasing circumferential angle. Figures B-1, B-2, and 5-3 
show a second pressure rise caused by the shock reflection at tie 


il 


| 
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fin. Figures B-4, Be5, and B-6 show a small initial pressure rise 
followed by a much larger rise. The small rise is caused by the 
iiffraction of the shock over the fin; while, the second one is 
caused by the influence of the reflected shock on the flow. The 
fin can be expected to have the following influence: 


l. It will cause the pressute on the fin side to be 
larger than that on the finless side. 


It will locally retard the shock on the fin side. 
The retarding of the shock will cause the finless side to be 
engulfed sooner. 


III. DESCRIPTION OF PRESSURE FIT 


In order to integrate the pressure on the surface cof the 
missile assembly at a given time, +, it is desirable to express 
the pressure as a function of coordinates “f the assembly surface 
at any desired time. Two convenient coordinates for this purpose 
are shown in Figure 7. Coordinate s is the longitudinal distance 
as measured along the surface from the arsembly tip, and & is the 
circumferential angle as measured trom the windward side. Once 
the pressure as a function of these coordinates is found, the 
forces on each section can be calculated from the integral, 


F= sf p(s,)aa, (2) 
A(s,6) 


where: 
F is the force on a given section, 
A(s,@) ie the sectional area engulfed by the blast wave. 


The empirical function, p(e,9), should be restricted to one 
which is periodic in @ with a period of 27 radians. It is also 
desirable that the function be easily integrable, so that 
numerical methods jo not have to be resorted to. If symmetry is 
assumed, then the function mist be even in 9. Thus, only one aide 
of assembly has to be considered, i.e., OS 9 <7. A cosine series 
in @ has these desired properties. Expressed as a cosine series, 
p(s,@) has the form 


L 
(8,6) = £  By(e) cos [(3-1) a), (3) 
j= 
where B, (s) are coefficients. In order to keep p(s,6) simple, 
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B:(s) are taken as polynomials ins. Thus, 


8 


By(s) bet? . (4) 
Phe complete function is 
K L 
p(s,6@) = ae ae cos [ /J-1) a]. (5) 
tel j=l 


Many other functions might be chosen to describe the preesure 
distribution. However, the one above has all of the desired 
properties, i.e., easily integrable, even in 9, and periodic with 
a@ period of 27. The coefficients by are determined by fitting 
the pressure data to the function at a given time, t, by the 
method of least squares. The summation limite, K and L, depend 
upon the number of data points available, i.e., the product of K 
and L must not exceed the number of data points available at 
time, t. 


The total caiculation time for each section is, for conveni- 
ence, divided into tw. periods, the engulfment period, and the 
post-erculfment period. The engulfment period is the time 
requir:d for the shock to engulf a given section. The post- 
engulfment period is the time after engulfment. During the 
engulfment period, the number of available data points varies 
because the number of transducers engulfed by the shock wave 
changes. Thus, during this phase, the limits, K and L, change 
from time to time. During the post-engulfment phase, the limits 
remain constant at K = 3 and L = 6. Of the many combinations 
considered, these values of K and L were found to give the best 
results. 


The calculation times during the engulfment period must be 
chosen with care. The reason for this is shown with the aid of 
Figure 8. Consider the diffraction of the shock around a parti- 
cular section of the missile assembly with the pressure trans- 
ducers located every 30° from the windward side. When the shock 
is between two transducer locations, the function p(s,9) may not 
describe the pressure accurately in the region between the trans- 
ducer last engulfed and the shock front. For example, if the 
shock were at 117°, only four data points would be available for 
fitting the pressure to p(s,6). Thus, there would be very little 
confidence in the function for angles greater than 90°. To avoid 
these inaccuracies, calculation times were chosen, during the 
engulfment period, as close as possible to the shock times of 
arrival at the various gauge locations. 
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IV. DESCRIPTION OF FORCE CALCULATION 


With the surface pressure given at a particular time as a 
function of s and 9, the force on a given section can be calculated 
by the expression, 


F=2 ff p(s,@)aa. (6) 
A/2 


During the engulfment period, A is the area of the section that 
has been engulfed by the shock. During the post engulfment period, 
A is the total area of the section. During the engulfment period, 
A is bounded vy the curve formed by the intersection of the shock 
front with the surface of the section. Of the three m’ -sile 
assembly sections of interest, two are conical frustrums (sect:ions 
2 and 3) and one is a cylinder (section 4). The intersection of 
the shock front with these two different types of sections is 
shown in Figure 9. 


For a cylindrical section, the equation for the curve of 
intersection is 


Yg = Te cos @ (7) 


where, Y, is the perpendicular distance of the shock front from 
the axis of the cylinder, r, is the radius of the section and 6 is 
the circumferential angle at which the shock front intersects the 
cylinder. Y, is positive when the shock is engulfing the front of 
the cylinder, but changes sign after the shock crosses the center 
of the cylinder. The corresponding relation for the frustrum 
rection is 


- —# 5 
Ys = Sing °° 6 (8) 


where R is the length of a ray from the theoretical apex of the 
frustrum to a point on the boundary cur'’e, and a is the half angle 
of the extended cone (see Figure 9b). 


Before the force calculations are described, three different 
gets of functions, f14, @1j, and hy4, will be defined. These 
functions occur naturally as the results of recurring integrals, 
and represent the (ij) terms in double summation expressions. 

The first function along with the defining integral is: 
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(9) 


(9a) 


(9b) 


(10) 


(10a) 


(10) 


(10c ) 
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(10e) 
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The final function ise: 


Y2 Xi ie 
hyy(yi,y2,x) = f Sf (x-d)*"" cos ( (1-2) y ] xdxay (11) 
Yl Asec y 
or: 
oe 
bap = pg? Coe) xa -aleyyr 1) (11a) 
n=0 
to” Cin [ vole ) ~ ab ) 
hyy = y — (= sin Jy - A’&rg j>.l (1lb 
7 n-0 1 J 4 


In the equations for hy}, 
I= i-n+l 


Je jel 


Cin = (i-2) a" / (d-1-n) to 
A= Y,g/sin a 


The parameter, d, is a reference distance whose physical significance 
will be explained later. 


Since the force calcuiations for the frustrum sections are more 
complicated than those for the cylindrical section, thé two cares will 
ve treated separately: 


A. Cylindrical Sections 
1. Engulfment Period 


The force on the cylindrical section at some time, t, 
during the engulfment period is the integral of the pressure function 
over the area bounded by the curves: 


@ = are cos (Ye 7he)) (12) 

8 = 69, (12a) 

; 8 = 8], (12d) 
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are the longitudinal distances of the end of the 
igure 10). Thus, the force is calculated by: 


85 5 
F= s¢ +s rep(s,@)d@ds « 
S$] @) 
Substituting equation (5) for p(s,9) into equation (23), 


K G 50 o 4 
F =r } Eb ST si-l cos C(-1)e) 
i=l j=l $3, 0 
xr, using equation (9), 
K L > 
F = 2r, Nee fe5n)4 (sj, So, 9) (1ha) 


2. Post Engulfment Period 


During the post engulfment period the pressure function is 
integrated over the total sectional half area. Since for this case, 


a 
@ « 
- , 


rx 
ae) 
"3 


by3 fy-1,5 (81> 82» 7) (15) 


B. Frustrum Sections 


a a Engulfment Period 


For a frustrum section, the area swept out by the shock is 
more easily visualized if the section is folded out. Since only the 
sectional half area is being considered (symmetry assumption), this 
area can be represented by 2 portion of a sector whose angle is 
B. = ™ sin a (see Figure lla). 


The shock boundary curve is then represented as a curve 
on the surface of the sector. The equation of this curve is, by 
equation (3), 


R = A sec 6 =Asec (f/sina ) 


The second equality gives the boundary curve in terms of the polar 
coordinates R and ®; where, 6 is defined as 6 sin a. A folded out 
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frustrum section is shown in Figure 1l for various stages of engulf- 
ment. The subscript, 1, refers to the small end of the section; while, 
the subscript, 2, refers to the large end. The quantities, 

r, = R; sin a and ro = Ro sin q , are the radii of the small and large 
end of the: section, respectively. The quantities: 


Oy = are cos (Yg/r)) , (17) 
4, = are cos (Y,/ra) , (18) 
Ru = A sec OL. (19) 


define the end points of the shock boundary curve. Equations (17) and 
(18) can be written in terms of the polar angles, Fy and Py, by multi- 
plying these equations by sin «a. The polar coordinates of the end 
points of the boundary curve are given in Figure 11 for each stage of 
engulfment shown. Since the pressure function, p(s,9) is dependent 
upon s, and since the area of a frustrum section is most conveniently 
expressed in terms of R, the relation between s and R must be used to 
express both quantities in terms of the same variable. This relation 
is, 


Rews4+a4, (20) 


where d is defined as, R) - 8}. 


The force calculations for each of the engulfment stages 
of Figure 11 are described below: 


a. Po< %Yg eT] 


Referring to Figure lla, 


By Ro 
r=2/ / p(s,6) Rade (21) 
O R 
By Ro 
F=2 sina /{ f p(s,6) RdRdO. (21a) 
O R 


d). -r, < yg<o @).-1 < ¥g< —F, 
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Substituting equation (5) for p(s,@) gives, 


kyu ; Ro er 
F = 2 sina tay ff, (Rea)? cos [(j-1)9) rarae (210) 
i=] jel O Asec 6 
or using equation (11), 


| ome 
FS Qatar 2%: yy (0, e ugha) (21c) 
fel jel 


Bee Lg ery 


Referring to Figure 1lb, 


one 


61, Ro Py Ro 
F=2sina}ys / p(s,@)RaRa9 +s f _ p(s,@)RdRAA (22) 
1@) Rj On, Asec 9 


Using equations (5) and (20), gives, 


K L u Ro r 
F=2sina¢ 5 £F byfl f si) cas ((s-1)8 ] (s+d)dsd9 (22a) 
f=1 j=l O Ry 
LS 
Yu Re re \ 
+ 7,  _(Rea)*™ cos { (j-1)0] RaRaa, >. 
Gy, Asec 6 | 


The first integral can be written in terms of fy; by equation (9). 
The second integral can be written in terms of hj; by equation (10). 
Thus, 


K L z 


F = 2 sina ne a by3 | £15 (GL» Rr, Ro) + dfy.1, 5 (Oz, Ry» Ro) (22b) 


Ce re) | 


2k 


ce. Yg = O 


Referring to Figure lle, 


1/2 Ro 
F =? sin. ys fs p(s,9) RAaRdO (23) 
fe) Ry 


or using equations (5), (20) and (9), 
K L 


2 sin diy (Py) (1/2s Ry» Ro) + Af yy, y(n/2, Ry, Ro) J) (23a) 
t=1 j=l 


es ee 


Referring to Figure lld, 


Dy Ro 9 Ry 7 
F=2sinalf f p(s,@)RaRa9 +f Sf _ p(s,6)RaRdo |. (2h) 
9 Ry Oy Asec 9 


Using equations (5), (20), (9), and (10) gives, 


K L 
F=ee2sin@ ¢ a by (115 @,, Ry» Ro) + dfs. (Oy, Rj, Ro) (2ha) 


im ‘= 


el je 
+y5 (1, Fy, Ri) 
tS ea aac | 
Referring to Figure lle, one sees that Py = ™/sin a 


and that 41 = 7, thus, 


‘ 


K L 
F = 2 sin & Qj by (P15 (Ou» Rr» Ro) + dfy-7,5(8y, Ris Ro) (25) 


+hay(7, By, Pr) ] 


3Recall that Y, changes sign after the shock passes the center of the 
section. This is true by the definition of ¥,. 
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. Host Engulfment Period 


During the post engulfment period, the pressure intcgra- 
tion {s performed around the total half sectional area. Thue, the 
force 18 calculated by the equation, 


n Ro 
2 sina ! s p(s,@) RdRdu (26) 
9 Ry 


yr using equations (5), (20), and (9), 


K L < 
Fe 2ein4 by 3 (tis) Ry, Ro) + dfy.105(", Ry, Ro) } (26a) 
{<1 je] 


Therefore, by using equations (21) through (26), the forces 
for a given frustrum section can be calculated for both the engulfment 
and the post engulfinent periods. 


C. Lescription of Computer Program 


The above results were used {in a computer program for calcu- 
lating the forces on given sections of the SPARTAN missile assembly. 
[his program, named PREDIN, is written in extended FORTRAN IV for the 
NWL CDC 6700 computer. A basic flow chart of this program is given in 
Figure 12. [his program takes pressure vs time data at given trans- 
ducer locations in-the form of punched cards, determines the pressure 
at each location for a given time, reads and appropriate limits K and 
L for the double summation of the pressure function, fits the data to 
the pressure function, integrates the obtained function over the 
appropriate sectional areas, and finally plots force-time histories for 
each missile assembly section of interest. 


V. RESULTS 
A. Results of Calculations 


Using PREDIN, force calculations were obtained on sections 2, 
3 and 4 of the missile assembly (see Figure 1) for DASACON tests 114- 
118. The calculations were performed on each section, at given times, 
for a total period of about seven milliseconds. Table 2 shows the 
calculation times for the engulfment period and the cc:responding 
values of K and L for cach test. Calculations were performed at 
intervals of .1 milliseconds during the post engulfment periods. . 
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READ PRESSURE ~- TIME DATA 


——— ee 


DETERMINE PRESSURE AT SPECIFIC 
LOCATIONS FOR A GIVEN TIME 


FIT PRESSURE FOR A GIVEN TIME 
TO FUNCTION, P(SO) 


INTEGRATE FUNCTION OVER GIVEN REGIONS 


CALCULATION 
PERIOD 
1S 
COMPLETE 


NO 


YES 


PLOT FORCE VS. TIME 


FLOW CHART FOR COMPUTER PROGRAM, PREDIN 
FIGURE 12 


mea YS 


Times 
Teast 14 


(msec) 
.150 
.550 
. 300 


Times 
“est 15 


_(mnec) 
50 


"Imes 
Pest 16 


(msec) 
-160 


TABLE 2 
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Times 
Test 17 
msec 


VALUES OF SUMMATION LIMITS FOR VARIOUS CALCULATION TIMES 


Times 

Test 18 
(msec) 6 L 
-200 ] 2 
-600 “ 2 
850 e 2 
1.000 2 2 
1.150 3 3 
1.750 3 hk 
2.350 3 ? 
2.950 3 6 


Ail times are referenced to the shock time of arrival at station Py-0 
see Figure 4). Tables of all of the byy coefficients of p(s,@) are 
given for each calculation time of each test in Appendix C. Appendix D 

shows comparisons between pressures calcuiated by p(s,9) and experi- 
mental data. ‘The appendix contains experimental data from test 114 and 
16. The comparisons for }14 are typical of tests 114, 117 and 113, 
and the comparisons for test 116 are typical of tests 115 and 116. 
Four times were choren for these comparisons, two during the engulfment 
period and two during the post engulfiment period. The solid lines in 
the figures were calculated by p(s,9). The points are experimental 
data. Graphs of force vs time on sections 2, 3, and 4 of the missile 
assembly are given for each test in Appendix E. These curves are 
characterized by an initial fast rise, followed Ly a nore gradu] rise 
a peak value, followed by a very gradual decay. The peak value 
‘curs 4t a time approximately equal to the engulfment period of the 
section. 


Effect of Fin 


As discussed enrlier, these results were obtained by assuming 
& symmetric pressure distribution. The effect of the fin on the 
pressure distribution of section 3 (control section) can be shown from 
the data of tests 115 and 116. These tests were conducted at approxi- 
mately the same incident shock pressure (see Table 1). Because of the 
rotation of the missile assembly after test 115, the two tests gave 
pressure measurements at 30°, 60°, 75°, 105° and 150° on the fin side 
of section 3. Mgures 13 through 16 show pressure data from the fin 
side vs angle compared with corresponding data on the finless side. 
The comparisons were made at times of 1.50, 2.00, 2.50 and 4.00 milli- 
seconds, respectively. The data from the fin side are connected by 
dashed lines; while solid lines are draw through the averages of the 
data from the finless side. These lines are drawn only to facilitate 
data reading. These figures clearly show the increased pressure due to 
the effects of the fin. These effects are smaller for section 2 and 4 
which do not have a fin. An estimate of the percent differences 
between the forces calculated using the symmetry assumption and that 
of the actual pressure distribution can be obtained by the following 
procedure: 


The force on section 3, at a given time, is estimated by 
multiplying the average of the pressure at the given time by the 


appropriate area. Thus, the force estimated without aceuming a 
symmetrical pressure distribution is 


Fy (p+ B) Afe ; (27) 
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while, the force estimated with the symmetry assunption is 


F = pA (28) 


where, Pp is the average pressure on the fin side, and Pp is the 
average pressure on th: finlees side. A is the: <otal area swept out 
by the shock wave. The ratio F,/F is given by, 


F,/F = (pp + B)/ep = (1 + By/p)/2 (29) 


This ratio is plotted against time in Figure 17. The maximum vuilue of 

the ratio occurs near the engulfment time. Figure 17, therefor:, gives 
as estimate of the effect of assuming a symmetrical pressure distribu- 

tion in the force calculations for section 3. This effect is expected 

to be less for sections 2 and 4 since they do not contain a fin. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


The results of this report show that the pressure distribution at 
a given time on the surface of the SPARTAN missile assembly can be 
represented by the function, 


aa 
p(s,0) = 5 5 by3 st") cos [(j-1) @ J 
i=1 jel 


The coefficients and the number of terms in this representation depend 
upon the amount of experimental dita available. For determining the 
forces on the missile assembly, the number of data points obtained 
during the DASACON tests were adequate. Analyses which require a more 
accurate definition of the pressure distributicn would, in general, 
require more data. 


The term, "force", used in this report applies to the integration 
of the above equation over the appropriate area. The actual directional 
forces on the various assembly sections could be obtained with the data 
presented and described in this report. For example, the directional 
forces on a given assembly section can be determined by, 


Fa = sf s -p(s,6)adaA, (30) 
A(s,@) 


where n is a unit vector normal to the surface. Vector n, can be 


expressed in terms of an orthogonal set of vectors A, N and K, shown in 
Figure 18 for a frustrum section. Unit vector A is directed along the 
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UNIT NORMAL ON THE SURFACE 
OFA 
FUSTRUM_ SECTION 


figure |8 
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sectional axis, N is directed perpendicular to the axis in the direc- 
tion of the blast wave propagation, and K is directed out of the plane 
formed by A and N. Thus, n can be written as, 


n= cos 0 A - sing cos # N + sin a sin 9 K Px) 


Since p(s,9) 1s symmetrical about the plane perpendicular to K, the 
total] force in the K direction must be zero. Substituting equa- 


tions (5) and (31) into equation (30) gives: 

ee mil : 

F, = sin « » byy ys 88) cos [(J-1)@]cos 9 4A (32) 
fet jek A(s,9 
K aC 

Fi cog o" 3 yr, by) Ss git) cos ((3-1)9] dA (33) 
$-l j=l A(s,9) 

Fy = 0 (34) 


"hese are the vector force components for a frustrum section in the 
axial, normal and K directions respectively. Similar expressions for 
a cylindrical section are obtained by setting « = O in the above 
equations. The same integration Jimits for the various stages of 
engulfment and the post engulfment period given in section IV can be 
used in evaluating the above integrals. Solutions to these integrals 
can be found in standard tables, although care must be takua t> insure 
that all possible cases are considered. 


The test results indicate that the presence of a fin on one side 
of the control section may cause the force on that section to be as 
much as 1.3 times as great as that encountered without a fin. ‘This 
result indicates that a symmetrical distribution should not be assumed 
when a high degree of accuracy in the force calculations is desired. 
The same equation for fitting the data can be used for both the 
symmetrical and the asymmetrical distribution. However, the fit for 
the asymmetrical distribution requires more data to be taken on the 
fin side of the assembly. For this reason, it is recommended that, in 
the future, an equal number of pressure gauges be placed on each side 
of asymmetrically loaded test items. 


NWL possesses a large amount of surface pressure data obtained 
from various blast loading tests on sectional assemblies of the Sprint 
and tthe Nike Hercules missiles, in addition to the SPARTAN data dis- 
cussed in this report. These missile assemblies vary in size, but have 
the same basic shapes, i.e., cylinders, frustrums of cones. These data 


ot 


include the results of tests performed over a wide range of incident 
shock pressures and missile orientations. It is recommended that a 
thorough aerodynamic and statistical analysis of the data be made. 
This analysis should lead to the formulation of an empirical pressure 
distribution model for the air blast loading of objects similar in 
shape to the missile assemblies. 
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SMOOTHED PRESSURE-TIME HISTORIES 
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